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Abstract 

Rotational and rovibrational spectra of methyl cyanide were recorded to analyze interactions in low-lying vibrational states and 
to construct line lists for radio astronomical observations in space as well as for infrared spectroscopic investigations of planetary 
atmospheres. The rotational spectra cover large portions of the 36-1627 GHz region. In the infrared (IR), a spectrum was recorded 
for this study in the region of 2vg around 717 cm -1 with assignments covering 684-765 cm -1 . Additional spectra in the vg region 
were used to validate the analysis. 

Information on the K level structure of CH3CN is almost exclusively obtained from IR spectra, as are basics of the J level 
structure. The large amount and the high accuracy of the rotational data improves knowledge of the J level structure considerably. 
Moreover, since these data extend to much higher J and K quantum numbers, they allowed us to investigate for the first time in 
depth local interactions between these states which occur at high K values. In particular, we have detected several interactions 
between vg — 1 and 2. Notably, there is a strong Avg = ±1, A K = 0, AI = +3 Fermi resonance between vg = 1 1 and vg = 2 +1 
at K — 14. Pronounced effects in the spectrum are also caused by resonant Avg = ±1, A K = +2, AI = +1 interactions between 
us = 1 and 2 at K = 13, l = -l/K = 11, l = 0 and at K = 15, l = +1/K =13,1= +2. An equivalent resonant interaction occurs 
between K = 14 of the ground vibrational state and K = 12, / = +1 of vg = 1 for which we present the first detailed account. A 
preliminary account was given in an earlier study on the ground vibrational state. Similar resonances were found for CH3CCH and, 
more recently, for CH3NC, warranting comparison of the results. From data pertaining to vg = 2, we also investigated rotational 
interactions with 14 = 1 as well as Avg = ±1, AK = 0, AI = +3 Fermi interactions between vg =2 and 3. 

We have derived N 2 - and self-broadening coefficients for the vg, 2vg - v x , and 2vg bands from previously determined V 4 values. 
Subsequently, we determined transition moments and intensities for the three IR bands. 
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1. Introduction 

Among the more than 180 different molecules detected in 
the interstellar medium (ISM) or in circumstellar envelopes 
(CSEs) of late type stars, there are 56 with a CN moiety, see, 
e.g., the Molecules in Space pageQof the Cologne Database for 
Molecular Spectroscopy, CDMS mm or the Bibliography of 
Astromolecule^ of the Astrochymist web-page]^] Almost two 
thirds of these can be viewed as cyanides in a wider sense. 
The alkyl cyanides n-propyl cyanide, / 7 -C 3 H 7 CN, also known 
as 1-cyanopropane or butyronitrile, and Ao-propyl cyanide, i- 
C 3 H 7 CN, also known as 2-cyanopropane or Ao-butyronitrile, 
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both recently detected in the very massive and luminous galac- 
tic center source Sagittarius B2(N) (Sgr B2(N) for short) [3|[4| 
are the largest complex organic molecules detected in the ISM 
thus far, and, after HCnN, the second largest cyanides. More¬ 
over, Ao-propyl cyanide was the first branched alkyl compound 
detected in space |[4]]. Ethyl cyanide, C2H5CN, the next smaller 
homolog, was detected much earlier in the Orion Molecular 
Cloud as well as Sgr B2 0 . Its M C isotopologs were detected 
in Orion IRc2 | 6 ] and Sgr B2(N) |7). In addition, rotational 
transitions pertaining to the two lowest-lying vibrational states 
were detected in G327.3-0.6 | 8 ) and Sgr B2(N) (9). Transi¬ 
tions of higher excited states were detected in molecular clouds 
of Orion m and Sgr B2(N) IfTTl . 

Methyl cyanide, CH 3 CN, also known as acetonitrile or 
cyanomethane, is the smallest alkyl cyanide. It is a trace con¬ 
stituent in Earth’s atmosphere which was studied, e.g., by high- 
resolution spectroscopy using microwave limb sounding G2 
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and balloone-borne IR spectroscopy tm It is produced mainly 
by biomass burning mm. Acetonitrile is also a useful sol¬ 
vent in the chemical industry. 

Methyl cyanide is a very important molecule in space. First 
identified in Sgr A and B more than 40 years ago, it was 
among the first molecules to be observed by radio astron¬ 
omy m. It was also seen in comets, such as Kohoutek 
03, in the atmosphere of Saturn’s moon Titan fT 8 l . in dark 
clouds, such as TMC-1 on, around the low-mass protostar 
IRAS 16293-2422 lf20ll . in the CSE of the famous carbon-rich 
late-type star CW Leonis, also known as IRC +10216 ED, and 
in external galaxies, such as NGC 253 |22| . Several isotopologs 
of CH 3 CN were detected as well in the ISM, including the rare 
CFFDCN 1231 and possibly even ' ’CH^CN ITTl . Transitions 
of methyl cyanide within its lowest excited vibrational state 
i )8 = 1 were detected more than 30 years ago El- More re¬ 
cently, brief accounts on observations of transitions belonging 
to W 8 = 2 IfTTl [251 . (,'4 = 1 ITU as well as on transitions of 
vibrationally excited methyl cyanide with I3 C lUTl have been 
published. 

The identifications were based on extensive laboratory in¬ 
vestigations which began almost 70 years ago |26j. The most 
recent account on the ground state rotational spectrum of the 
main isotopolog CH3CN was given by some of the present au¬ 
thors along with analyses for several minor isotopologs ED; 
unlabeled atoms refer to 'H, 12 C, and 14 N. The CH3CN line 
list was based mainly on an accurate and extensive study l28l 
which employed earlier data for frequencies below 74 GHz 
ESlEIllIO- Information on the purely K -dependent terms 
cannot be determined by rotational spectroscopy for a strongly 
prolate symmetric top, such as methyl cyanide, unless from the 
very weak AK = 3 transitions that gain intensity from centrifu¬ 
gal distortion effects on the the dipole moment. Thus the purely 
/^-dependent terms A, Dk, and an estimate of IIk were deter¬ 
mined from AK = 3 ground state energy differences |[32l . 

The lowest excited vibrational state of CH3CN is the doubly 
degenerate CCN bending mode ug = 1 at a vibrational energy 
of 365.0 cnT 1 l33l . It was studied first in 1961 by rotational 
spectroscopy between 92 and 222 GHz ( J" = 4 to 11) Ii34l . 
Soon thereafter, more accurate data with partly resolved 14 N 
hyperfine structure between 18 and 148 GHz ( J" = 0 to 7) 
were reported ll35ll36S . Later, data were obtained near 370 GHz 
{]” = 19) E3 and, employing a laser-sideband spectrometer, 
in selected regions between 589 and 1395 GHz ( J” = 31 to 75) 
EH- Direct /-type transitions between the k x / = +1 compo¬ 
nents were also observed with 21 < / < 45 Il39ll40ll . 

To our knowledge, there was only one prior high-resolution 
IR study of the +8 band B33H , although there were earlier low- 
to moderate-resolution studies of vg or indirect determinations 
of the vibrational energy from higher vibrational bands. A res¬ 
onant interaction with A k = -2 and Al — +1 between v - 0, 
K — 14 and yg = 1 +1 , K = 12 at / = 43 was reported 
ED, but could only be analyzed approximately because only 
two slightly perturbed transitions in the ground vibrational state 
( J" - 38 and 39) were observed. 

The parallel 2+8 band at 716.7 cm -1 l42l was treated as an 
isolated state in that study even though a rather strong Fermi 


interaction of vg = 2° with u 4 = 1, the non-degenerate C-C 
stretching mode at 920.3 cm -1 (W 488 » 30 cm -1 ) was proposed 
from a force field analysis Il43l . This interaction has not been 
considered in high resolution analyses of vg = 2 thus far. Tran¬ 
sitions to vg = 2 ±2 became observable only through <722 interac¬ 
tion with vg = 2° Ii42l . A Fermi interaction between vi = 1 +1 , 
and vg = 2~ 2 at K = 12 and 13 was proposed to explain devi¬ 
ations in the high- K part of the +7 band Ii44ll . This interaction 
was confirmed in a later study mainly concerned with v 4 , +7 
and 3+g bands between 850 and 1150 cm 1 E3; in addition, a 
Ak = -1 and Al = +2 interaction between vg = 2 ~ 2 and y 4 = 1 
at K = 6 and 5, respectively was found for the first time along 
with several other perturbations among the higher lying states. 

The rotational spectrum of CH 3 CN in its vg — 2 state has 
also been studied between 18 and 112 GHz ( J" = 0 to 5) 
fflSSHZll, near 370 GHz ( J" = 19) ll37l and, employing a 
laser-sideband spectrometer, in selected regions between 787 
and 1273 GHz U" = 20 to 68 ) If38l. 

Despite the numerous studies of methyl cyanide in various 
spectral regions, there has been no detailed attempt to study its 
spectroscopy globally. If there were effects, which were caused 
by vibrational states outside of the polyad under investigation, 
they were only considered through the effective parameters, see, 
e.g.. Refs. EDSa. 

In the present work, we have recorded extensive rotational 
and rovibrational spectra in order to study vibration-rotation 
interactions not only within vibrational states or within vibra¬ 
tional polyads, but also between vibrational polyads. In com¬ 
parison to earlier data, we access higher J and K quantum num¬ 
bers, have a more complete coverage of the data and typically 
higher accuracies. The data permitted the first detailed analy¬ 
sis of a resonant interaction between v — 0 and vg - 1 and the 
first detections and analyses of interactions between vg = 1 and 
2 as well as between vg = 2 and 3. Preliminary results were 
described earlier Il48l . In this article, we present an expanded 
analysis of these vibrational states of methyl cyanide with i >8 
up to 2 and provide preliminary analyses of interactions with 
vibrational states pertaining to the next polyad. We calculate 
rotational and rovibrational spectra from our resulting analysis 
to support methyl cyanide observations in diverse environments 
in space as well as in Earth’s atmosphere and those of other 
planets or related objects. 

Interactions between vibrational polyads of symmetric top 
molecules, such as methyl cyanide, are facilitated by the com¬ 
paratively small vibrational energy of 365 cm ~ 1 for ug = 1 com¬ 
bined with the large Coriolis coupling between the two / com¬ 
ponents which causes levels with AK = 0, but different / to 
repel each other increasingly at higher K such that the lower 
energy / component can interact with the ground vibrational 
state at specific K whereas the higher energy / component in 
turn can interact with the different / components of i >8 = 2 also 
at specific K. Vibrational polyads are groups of one or more 
vibrational states separated by approximately equal amounts of 
vibrational energy, see e.g. Ref. Il49l for a recent review. 

Interactions in symmetric top molecules between vibrational 
polyads were studied first in the case of propyne, CH 3 CCH, also 
known as methylacetylene, which is isoelectronic to methyl 
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cyanide. A Ay = ±1, A K — 0, A/ = +3 Fermi interaction at 
K = 12 between tqo = U 1 and v\o = 2 +2 as well as a Coriolis 
interaction with A K + 1, A/ = +2 at K — 10 and 9, respectively, 
between tqo = 1 1 and vg = 1 +1 were first revealed mainly 
by rotational spectroscopy (50) and later improved by infrared 
as well as additional rotational spectroscopy (51) . The latter 
study also revealed a rotational (Coriolis-type) resonance be¬ 
tween vio = 2 and vs = 1 , which is the lowest vibration of the 
next higher tetrad. Interactions between the dyad and the tetrad 
were studied in greater detail later (52). More recently, such 
Fermi and a-resonances were analyzed between yg = 1 and 2 
in methyl isocyanide l53l . The vibrational energy of = 1 is 
lower, only 267.3 cm -1 , such that the Fermi interaction occurs 
already at K — 9 and 10, while they occur at K — 12 and 14 for 
the corresponding resonances in methylacetylene and methyl 
cyanide. 

We are not aware of resonant interactions between the ground 
and an excited vibrational state in the case of a symmetric top 
molecule other than the case of CH 3 CN. However, a Coriolis- 
type interaction between the ground vibrational state of HNCO 
and its lowest v$ = 1 state was proposed and its strength esti¬ 
mated IS). An analysis was presented for the analogous case 
of the isoelectronic HN 3 molecule 15511 . A remarkable case is 
vinyl cyanide, C 2 H 3 CN, also known as acrylonitrile, for which 
the ground vibrational state was connected to its lowest vi 1 = 1 
state through A K = 5 and 6 interactions without observation 
of high-resolution infrared transitions (56). Recently, vu = 1 
was connected to 015 = 1 and vu = 2 as well as the latter 
two to each other (57). The density of vibrational states in¬ 
creases with energy, and therefore chances to connect states at 
lower quantum numbers and through interactions with smaller 
A K increase. However, observation may be difficult in practice 
because transitions in higher vibrational states are weaker than 
corresponding ones in lower vibrational states. In addition, the 
number of at least similarly strong rotational transitions per fre¬ 
quency unit increases also with vibrational state. 

2. Experimental details 

2.1. Rotational spectra at the Jet Propulsion Laboratory 

Most of the CH 3 CN rotational data were extracted from 
broad frequency scans taken with the JPL cascaded multi¬ 
plier spectrometer (58)■ Generally, a multiplier chain source 
is passed through a 1 - 2 meter pathlength flow cell and is de¬ 
tected by a silicon bolometer cooled to near 1.7 K. The cell 
is filled with a steady flow of reagent grade acetonitrile at 
room temperature, and the pressure and modulation are op¬ 
timized to enable good S/N ratios with narrow lineshapes. 
The S/N ratio was optimized for a highcr-K transition (e.g. 
K — 12) because of the very strong ground state transitions 
of the main isotopolog with lower K, which may exhibit satu¬ 
rated line profiles. This procedure enables better dynamic range 
for the extraction of line positions for rare isotopologs and 
highly excited vibrational satellites. The frequency ranges cov¬ 
ered were 440-540, 619-631, 638-648, 770-855, 875-930, 
967-1050, 1083-1093, 1100-1159, 1168-1198, 1576-1591, 


and 1614-1627 GHz, as in our previous analysis ED- In ad¬ 
dition, we obtained new data at 530-595 GHz and data with 
improved S/N in the 780-835 GHz region. Most of these 
multiplier sources were previously described (58). However, 
the multiplier chain with frequency range coverage between 
967-1050 GHz was not described in that work. This chain con¬ 
sists of two cascaded triplers after the amplified W-band stage, 
the peak output power is near 100 pW. The efficiency of fre¬ 
quency multipliers is usually strongly frequency-dependent. In 
addition, recording conditions and sensitivities of detectors can 
have strong influences on the quality of the spectra. Particu¬ 
larly good S/N ratios were reached around 600, 800, 900 and at 
1100-1200 GHz. The S/N ratios changed considerably within 
each scan and were usually lower towards the edges. Uncer¬ 
tainties of 50 kHz were assigned to isolated lines with good to 
moderate S/N ratios, smaller uncertainties, down to 10 kHz, in 
cases of very good S/N and very symmetric lineshapes, larger 
uncertainties if the S/N was poorer or if the lineshape was rather 
asymmetric, e.g. because of the proximity to another line. 

2.2. Rotational spectra at the Universitdt zu Koln 

All measurements at the Universitat zu Koln were recorded 
at room temperature in static mode employing different Pyrex 
glass cells having an inner diameter of 100 mm with pressures 
in the range of 0.1-0.5 Pa at lower frequencies, increased to 
around 1.0 Pa at higher frequencies. The window material was 
Teflon at lower frequencies, whereas high-density polyethylene 
was used at higher frequencies. Frequency modulation was 
used throughout with demodulation at 2 /, causing an isolated 
line to appear close to a second derivative of a Gaussian. 

The J - 2-1 and 3-2 transitions of several vibrational states 
around 37 and 55 GHz, respectively, were recorded with an 
Agilent E8257D microwave synthesizer as source and a home- 
built Schottky diode detector. A 7 m long double pass ab¬ 
sorption cell was used for these measurements. All other mea¬ 
surements were carried out in shorter (around 3 m) single pass 
cells. The / = 4 - 3 transitions of ug =2 around 74.0 GHz 
were recorded using a backward-wave oscillator (BWO) based 
4 mm synthesizer AM-MSP 1 (Analytik & Mefltechnik GmbH, 
Chemnitz, Germany) as source and a Schottky-diode as detec¬ 
tor. A small number of methyl cyanide rotational transitions 
were recorded around 850 GHz and around 880 GHz with the 
Cologne Terahertz Spectrometer ll59l using a BWO as source 
and a liquid helium cooled InSb hot-electron bolometer (QMC) 
as detector. Further, extensive measurements were carried out 
in parts of the 1330-1501 GHz region employing a frequency 
multiplier (Virginia Diode, Inc.) driven by an Agilent E8257D 
microwave synthesizer as source and the InSb bolometer as de¬ 
tector. This spectrometer was used previously for investigations 
of CH 3 SH and 13 CH 3 OH (60]|6D- 

2.3. Infrared spectrum 

One infrared spectrum of CH 3 CN was recorded between 600 
and 989 cm -1 at 0.0016 cm -1 resolution using an MCT detec¬ 
tor with a Bruker 120 HR Fourier transform spectrometer at Pa¬ 
cific Northwest National Laboratory (PNNL); the 2vg region is 
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Figure 1: The 2v% spectrum of CH 3 CN recorded at 0.0016 cm -1 resolution 
using a Bruker 120 HR FTS at PNNL; the optical path is 19.2 m, and the gas 
pressure is 0.226 Torr at 293 K. The sample included ~0.4 % CO 2 and ~3 % 
OCS for frequency calibration. 


731 731.2 731.4 731.6 731.8 732 



Figure 2: Comparison of CH 3 CN lines at 731.5 cm -1 retrieved by least squares 
curve-fitting. The lower panel has the observed spectrum in black overlaid by 
the synthetic spectrum in red. The upper panel shows the residual difference 
(calc.-obs.) in percent. The optical path is 19.2 m, and the gas pressure is 
0.226 Torr at 293 K.The vertical lines between the two panels mark the mea¬ 
sured features (every 10 th line is taller). 



Figure 3: Model of the methyl cyanide molecule. The C atoms are shown in 
gray, the H atoms in light gray, and the N atom is shown in blue. The a-axis is 
along the CCN atoms and is also the symmetry axis. 


shown in Fig. [Tj For this, a multi-pass absorption cell, set to an 
optical path length of 19.20 m, was filled to 0.226 Torr (30.1 Pa) 
of CH3CN at 293.0 K. Small amounts of OCS and CO2 were 
added to the sample for frequency calibration; comparison of 63 
well-isolated lines of the V 3 fundamental of OCS at 860 cm 1 
f62l produced a calibration factor 1.000000868 ( 21 ) with an 
rms of 0.0000182 cm -1 . Nearly 12000 line positions and rela¬ 
tive intensities between 690 and 765 cm -1 were retrieved using 
non-linear least-squares curve-fitting 11631 . However, the corre¬ 
sponding precisions in the 2 vs region were an order of magni¬ 
tude worse than the OCS calibration because the observed spec¬ 
trum was so congested. Fig. [2] shows the quality of the curve 
fitting for a small section of 2 vg. 

Two other room-temperature spectra spanned the vg region 
(250 to 400 cm -1 ) at 294 K using pressures of 0.78 Torr 
(103.9 Pa) with a 19.2 m path and 1.01 Torr (134.5 Pa) with 
a 3.24 m path. These were calibrated using residual water in¬ 
side the FTS l64l to obtain a calibration factor of 1.000000207 
(495) with an rms of only 0.000140 cm -1 . The spectra were 
very congested and no attempt was made to retrieve postions 
and intensities by curve-fitting. Previously reported CH3CN 
line centers If33l had better precisions in this second region, and 
so these new vs data were used mainly to search for hot bands 
and to confirm assignments in the 2v$ region. 


3. Results 

3.1. Overview of CH\CN spectroscopy 

As shown in Fig. [3j the light H atoms are the only ones not 
on the symmetry axis of CH3CN. Hence, methyl cyanide is 
a strongly prolate molecule (A » B). Rotational transitions 
with A K = 0 are strong because of the large dipole moment 
of 3.92197 (13) D ll65l . A K = 3 transitions only gain inten¬ 
sity through centrifugal distortion effects and are usually very 
weak. In the absence of perturbations, the purely axial param¬ 
eters A (or A - B), Dk , etc. cannot be determined by rota¬ 
tional spectroscopy. Moreover, even rovibrational spectroscopy 
yields, strictly speaking, the differences A A (or A (A-B)), AD k, 
etc. from single state analyses, but not the ground state param¬ 
eters. In the case of CH3CN, analyses of three IR bands involv¬ 
ing two doubly degenerate vibrational modes, vg, V 7 + V 8 , and 
V 7 + vs - vg, permitted the axial parameters to be determined 
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Table 1: Energies (cm 1 ) and symmetries Sym of low-lying vibrational states 
of methyl cyanide. 


E / cm 1 
a 

1600 - 


State 

Sym. 

Energy 

Reference 


u = 0 

A 

0.0 

per definitionem 


Vg = 1 

E 

365.024 

this work 

1200 - 

vg = 2 ° 

A 

716.750 

this work 


vg = 2 2 

E 

739.148 

this work 

1000 - 

U 4 = 1 

A 

920.290 

Ref. |45i 


Vi = 1 

E 

1041.855 

Ref. (451 

800 - 

vg = 3 1 

E 

1077.786 

Ref. |45j 


vg = 3 3 

A 

-1122.15 

Ref. |45l 

600 - 

V4 = vg = 1 

E 

1290.0 

Ref. |69l 


V3 = 1 

A 

1385.2 

Ref. (69) 

400 - 

V~J — Vg — 1° 

A 

1401.7 

Ref. (69) 


V~! — Vg — l 2 

E 

1408.9 

Ref. (69) 

200 - 

Vg = 4° 

A 

-1426. 

this work 


< 4f- 

11 

00 

Si 

E 

1447.9 

Ref. (59) 

0 -1 


Vg = 1 

E 

1449.7 

Ref. (69) 

Vg = 4 4 

E 

-1514. 

this work 


15 


14 -- 14, 

tr i i ■ ; 


< ->- 12 



10 


11 

10 


5 


0 


5 



0 


v= 0 


/ = +1 / = -1 
Vo=\ 


/ = + 2 


/ = 0 
v r = 2 


/ = - 2 


10 
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Methyl cyanide has four totally symmetric and four doubly 
degenerate vibrational modes. The three lowest energy funda¬ 
mentals are the doubly degenerate CCN bending mode vg at 
365.024 cm -1 l33lk the totally symmetric CC stretching mode 
v 4 at 920.290 cnr 1 03], and the doubly degenerate CH 3 rock¬ 
ing mode v 7 at 1041.855 cm 1 03. They are comparatively 
weak with integrated room temperature cross sections in the 
range 2 to 5 x 10 19 cm/molecule ||66l [671. Interestingly, the 
2 v s overtone band is slightly stronger than v 8 or V4 06] [67]. 
In the case of V4, line-by-line intensity measurements have also 
been published 155] . A list of low-lying vibrational states, up 
to ~1500 cm -1 , is given in Table [l] 

The fairly large value of A , ~5.27 cm -1 , leads to a rapid 
increase in rotational energy with K which is shown up to 
~1700 cm 1 in Fig. [ 4 ] The highest K levels observed involve 
K = 21, and the J = K = 21 energy is at 2313 cm -1 . 

The vibrational energy of the lowest, doubly degenerate vi¬ 
brational mode y 8 = 1 is at 365.024 cm -1 1331 , similar to the 
J = K = 8 energy in v - 0 (Fig.[4|i. The Coriolis interaction be¬ 
tween the two / components of u 8 = 1, described by £ » 0.877, 
is large, close to the limiting case £ = 1. It effectively shifts the 
levels of the / = +1 ladder to lower energies and the energies 
of the I — -1 ladder to higher energies. As a result, K - 1 of 
/ = +1 is actually lower than K - 0. More importantly, levels 
with A K — A/ = 2 are near-degenerate, leading to pronounced 
<722 interaction within u 8 = 1. The interaction is highlighted in 
Fig.^for K = 4 of vg = 1 +1 and K = 2 of vg = 1 _1 . As a con¬ 
sequence of the Coriolis interaction within vg = 1, the K level 
spacing in ug = 1 +1 is considerably smaller than in v - 0 while 
it is larger in y 8 = l -1 . Simeckova et al. PTFI found evidence of 
a resonant Ai > 8 = ±1, AK = +2, Al — +1 interaction between 
v = 0, K = 14 and ug = 1 +1 , K = 12 (Fig. [5} at J = 43. They an¬ 
alyzed the interaction approximately because only two slightly 
perturbed transitions in the ground vibrational state (7" = 38 
and 39) were observed. 

Other global interactions within an isolated degenerate vi¬ 
brational mode of a strongly prolate rotor are rare, in contrast 


Figure 4: K level structure of methyl cyanide displaying the J = K energies for 
v% = 0, 1, and 2 with their / substates. Levels with A symmetry are drawn with 
thick lines. Coriolis coupling within v& = 1 causes near-degeneracies of levels 
having A K = A/ = 2 and hence strong #22 interaction shown exemplarily for 
one pair of low K by a thin, dashed (blue) line. The same applies to i>8 = 2 ±2 , 
where a <722 resonance occurs between K = 4 and K = 2 of the / = 0 and -2 
components, respectively (thin blue left right arrow). The spacing between K 
levels is fairly large compared to the vibrational spacing causing a A K = 0 , 
A/ = 3 Fermi interaction at K = 14 between i»8 = l -1 and v% = 2 +2 , displayed 
by a thick (red) left right arrow. Also shown are A K = +2, A/ = ±1 interactions 
between i>8 = 0 and 1 as well as between v% = 1 and 2 (medium thick green left 
right arrow). 

to nearly spherical or oblate rotors, see, e.g. Ref. 1701 and ref¬ 
erences therein. 

The basic features of the K level structure of vg - 2 are sim¬ 
ilar to those of vg = 1 , the main difference being three / com¬ 
ponents, l = 0, ±2 instead of two. Pronounced t /22 interaction 
occurs also within ug = 2 , in particular between I = 0 and 
l = -2 at low K with a resonant interaction K = 4 of u 8 = 2° 
and K = 2 of vg = 2 ~ 2 , as shown in Fig. [ 4 ] In addition, the 
K level structure in vg = 2 +2 is irregular at low K with K = 2 
being lowest in energy and K — 3 still slightly below K - 0. 
Such irregular K level structure in states such as vg > 2 may 
lead to local perturbations within the vibrational state with, e.g., 
AK — +1, Al = +2 (q 2 \), AK = ±3, Al — 0 (e), etc. In the case 
of several different such local perturbations, indeterminacies, 
as investigated in Ref. ED and references therein, may need to 
be considered. A global effect may be caused by a AK = ±4, 
Al = ±4 interaction However, this effect is rather small 

in a strongly prolate rotor, and 4 was barely determined for 
uio = 2 of propyne ED. 

3.2. y = 0 

Recently, we have revisited the ground state rotational spec¬ 
trum of several isotopic species of methyl cyanide Il27l . Be¬ 
sides new data, ground state transition frequencies for the main 
species were mainly taken from Ref. Il28l with additional lower 
frequency data 129 30 .Tl |. The J and K quantum numbers 
extended to 89 and 21, respectively, see also Table [2] 
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Table 2: Maximum J and K values and number of rotational and IR lines in vibrational substates of methyl cyanide used in the present / previous study* 3 , with 
numbers retained in the present fit from previous data in parentheses and rms error for subsets of data. 



v = 0 

va = 1 +1 

Dg = I” 1 

^2° 

Vi = 2 +2 

yg = 2 2 

*Anax(rot) 

89/89 (89) 

88/75 (45) 

88/75 (8) 

88/68 (20) 

88/68 (20) 

88/68 (20) 

^max(rot) 

21/21 (21) 

19/13 (7) 

17/11 (7) 

19/12(12) 

13/8 (7) 

11/8 (6) 

no. of rot lines^ 

316/225 (196) 

596/c (58) 

527/c (37) 

424/85+5°' (18) 

546/73+6 rf (24) 

299/45+5 rf (9) 

rms error 6. 

0.859/0.911 

0.807/- 

0.796/- 

0.982/- 

0.915/- 

0.957/- 

/max(IR) 

-/- (-/ 

71/71 (71) 

74/74 (74) 

66/65 

53/37 

66/57 

^max(IR) 

7/7 (7) 

11/11 (11) 

11 /11(11) 

12/10 

2/4 

7/4 

no. of IR lines* 

5/5 (5/ 

836/g (836) 

861/g (861) 

935/944 (0) 

40/44* (0) 

197/122* (0) 

rms error 6, 

0.716/0.643 

0.826/- 

0.797/- 

0.771/- 

0.849/- 

0.727/- 


a Previous data from Refs. I27II33II42H and references therein for v = 0, v& = 1, and ug = 2, respectively. 

b Each blend of lines counted as one line. Hyperfine splitting in the present and previous data was considered in the present fit. Transitions with K = 0 of /-doubled 
states associated with ug = l -1 and ug = 2 +2 , respectively, in the present data set. Lines, which were weighted out, are not counted. Transitions between vibrational 
states were counted for higher vibrational state. 
c Only total number given as 187. 

d Numbers after the plus were given zero weight in Ref. 1421 . 

e No rms error was given in Refs. EH HD; we assume parameter uncertainties are based on standard errors, i.e. the rms error is 1.0 by definition. The rms errors 
for data from Ref. E3 were calculated from the experimental data available at http://www.astro.uni-koeln.de/site/vorhersagen/daten/CH3CN/CH3CN/gs/ 

J No individual A K = 3 loops were given in Ref. E2 We used the five A K = 3 splittings from Table II in that work with the reported uncertainties. 
s Only total number given as 1705. Slight difference to our numbers mainly because of blended lines, see subsubsection |3.3.2| 

* 52 assigned 2vs - vy hot band transitions are mentioned in Ref. 14221; see subsubsection|3 .4.3| 


The first objective and the main achievement of the present 
study with respect to the ground vibrational state was the de¬ 
tailed analysis of the resonance between v = 0 and vg - 1 
mentioned in the previous subsection. Absorption features per¬ 
taining to v — 0, K = 14, J" = 42 to 49 and to pg = 1 +1 , 
K = 12, J" = 41 to 45 and 47 to 49 could easily be assigned in 
our spectral recordings because of their intensities. Moreover, 
the perturbed transitions showed regular and mostly small dis¬ 
placements from predictions which did not take the interaction 
into account. The largest shifts were ~24 MHz and ~ 16 MHz 
for the / = 43 - 42 and 44 - 43 transitions, respectively. After 
these data had been included in the fit, two pairs of cross ladder 
transitions between the two vibrational states were predicted to 
have sufficient intensities for observation. The weaker pair with 
/ = 44 - 43 were identified, but the ones with J = 43 - 42, pre¬ 
dicted to be stronger by about a factor of 2.5 than the former 
ones, were blended with stronger lines. Fig. [5] shows a detail 
of the energy levels around the levels perturbed most by the in¬ 
teraction between v = 0 and pg = 1 as well as the observed 
transitions. 

The measurements in Koln between 1.33 and 1.5 THz closed 
a considerable part of the gap between ~ 1.2 and ~ 1 .6 THz in the 
previous data set Il27ll . We reanalyzed also the spectra between 
~1.1 and 1.2 THz taken at JPL and realized that the uncertain¬ 
ties were very similar to those of Ref. 1251 . albeit with a slightly 
better coverage and extending to slightly higher K. Therefore, 
we employed the JPL data in our present analysis. 

We also attempted to assign rotational transitions with A K = 
3. Their predicted uncertainties were small enough to iden¬ 
tify them easily if the distortion dipole moment was as large 
as 10 4 D. However, even then they could be blended with 
stronger methyl cyanide transitions, of which several are not 


yet assigned. No assignments could be made, leading us to 
conclude that the distortion dipole moment is probably smaller. 
It is also possible that the calculated uncertainties of these tran¬ 
sitions are too small. 


3.3. p 8 = 1 

3.3.1. Rotational data 

The rotational spectrum of CH 3 CN in its pg = 1 lowest vi¬ 
brational state {E v ib = 365.02 cm -1 ) was reported on several 
occasions ll35lE6ll37ll38l . The J" = 19 data Ii37l had residuals 
in part much larger than the reported 50-150 kHz uncertain¬ 
ties already in the initial report, which were confirmed here. 
Therefore, and because of their more limited K range, we omit¬ 
ted these data from our line list. The laser side-band measure¬ 
ments f38l were also quite limited in K and in accuracy (around 
1 MHz) and were superseded by our more accurate data. Part 
of the lower frequency data up to 150 GHz Il35l [36l was also 
superseded by our data, the remainder was retained in the fit 
with reported uncertainties. We also used in our fit direct /-type 
transitions between the k X l = +1 components 139] |40l . 

The interaction between u = 0 and pg = 1 was already de¬ 


scribed in subsection 3.2 Interactions of pg = 1 with any higher 
vibrational states have not been reported in the literature prior 
to the present investigation. 

Based on the studies on propyne mm, we expected an 
interaction between pg = 1 1 and p x = 2 +2 at K slightly higher 
than 12 because of the slightly larger vibrational energy dif¬ 
ference in methyl cyanide compared to propyne, 374.12 cm -1 
compared to 340.41 cm . Inspection of the energy levels of 
these CH 3 CN vibrational states indicated a resonant interaction 
to occur at K = 14, see Figs. |4] and [6] Using the value of 
the main Fermi parameter from the propyne study f50l[5ll was 
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V— 0 V g = 1 + 1 


Figure 5: Detail of the energy level structure of methyl cyanide around the 
perturbation between the ground state and ug = 1. Transitions observable in our 
spectral recordings are represented by anows; those with solid lines indicate 
transitions in the final line list, including two transitions between the vibrational 
states. Two further transitions between v = 0 and vg = 1 with significant 
intensities, but overlapped with stronger transitions are displayed by arrows 
with dotted lines. 



Figure 6: Reduced energy plot in the region of the perturbations between vg = 1 
and 2. A A«8 = ±1, A K = 0, A 1 = ±3 Fermi resonance occurs near 1450 cm 1 . 
In addition, two resonant Avg = ± 1, A K = +2, A/ = ±1 interactions are indi¬ 
cated near 1330 cm -1 . 



Figure 7: Fortrat diagram in the region of the resonant Avg = ±1, A K = +2, 
A/ = ±1 interactions between vg = 1 and 2. Observed transitions are indicated 
by larger symbols. 

a reasonable initial guess along with published parameters for 
vg = 2. Our final results suggest the interaction to be strongest 
at J = 91 and 92, however, the highest assignments presently 
reach / = 65 - 64 for either state at these rather high values of 
K. The magnitudes of the perturbations between vg = 1 1 and 
vg = 2 +2 showed strong variations with K, as expected, but also 
with /, most likely because the energy differences between the 
two vibrational ladders also change significantly with J. The 
perturbations exceed 1 MHz at / = 24 - 23, K = 10 and reach 
37 MHz at K — 14. At J = 65 - 64, perturbations in the spec¬ 
trum larger than 1 MHz were already seen at K — 5, and the 
shift exceeded 300 MHz at K - 14. 

The data set up to J — 28 - 27 revealed perturbations in 
K — 15 of vg = 1 +1 which were not accounted for by this 
Avg = +1, AK = 0, A/ = +3 Fermi interaction. They were 
caused by a Aps = ±1, AK = +2, A! — +1 interaction with 
K = 13 of vg = 2 +2 (Figs. @ and |6]), the same type of interac¬ 
tion as between v - 0 and vg — 1. It is shown in Fig. [7] that the 
strongest mixing occurs at / = 52 and 53. Among the four tran¬ 
sitions with J = 53 - 52, the iig - 1 transition was too close to 
the lower frequency limit of its experimental spectral window, 
the vg — 2 transition was observed with a total perturbation of 
1.8 GHz. In addition, both cross ladder transitions were ob¬ 
served, but one was too close to a stronger line to determine a 
reliable and sufficiently accurate transition frequency. In addi¬ 
tion, the / = 54-53 transition frequencies were determined for 
both vibrational states. 

The same type of interaction showed noticeable effects in the 
spectrum at higher J (up to 49) between K = 13 of ug = 1 1 and 
K = 11 of vg = 2°, see also Figs. gand[6] with the largest per¬ 
turbation at J - 60, as shown in Fig. [7] The regular rotational 
transitions were detected for both vibrational states in the case 
of / = 61 - 60, both cross-ladder transitions were detected for 
/ = 60 - 59. The total perturbations amounted to more than 
1.5 GHz for the former two transitions and more than 3.0 GHz 
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for the latter two. The remaining associated transitions were 
about an order of magnitude weaker and were not identified. 

We mention for completeness a Ai>s = ±2, AK = +2, A/ = +4 
interaction between K = 15 of t)8 = l _l and K = 13 of vg = 3 +3 
whose strongest effects are currently predicted at / = 73/74. 
Should the interaction cause any measureable perturbation in 
the spectrum, it should be very localized because of the high 
order of the interaction. A Avg = +2, AK = +1, Al = ±2 
interaction is possible between K = 26 of vg = 1 and K - 25 
of i>8 = 3 +1 , but these levels are so high in energy that they are 
unlikely to be accessed. 

The overall coverage of our «8 = 1 is indicated in Table [2] 
Besides the few cross-ladder transitions between vg — 1 and 2, 
the data are limited almost entirely to regular rotational transi¬ 
tions. Transitions between different l components may occur 
because of Al = AK = 2 mixing facilitated by q. Such transi¬ 
tions were predicted to be very weak, and only one transition 
was identified, so we view this identification as tentative. We 
also searched for AK = 3 transitions, but, as in the case of the 
ground vibrational state, were unable to identify any with suffi¬ 
cient certainty. 

3.3.2. The vs band 

An investigation of vs was reported in Ref. J33], and 
we use the line list from that work as indicated in subsec¬ 
tion |2.3| Lines with full weight were assigned uncertainties 
of 0.0002 cm -1 , except for the small number of ' Q\o lines, for 
which 0.0004 cm -1 were used because of their persistent larger 
residuals in the fits. Lines with weight 0.1 were assigned un¬ 
certainties of 0.0006 cm -1 . A few lines (around 10 to 20) with 
weight 0.1 were omitted either because of very large residuals 
or because the line was blended, but was the weaker line, and 
the stronger line was not in the line list. In other rare cases, 
blended lines (12 x 2 ) were treated as intensity-weighted aver¬ 
ages. The effects caused by these differences in the previous 
and present line list on the parameter values, uncertainties and 
on the quality of the fit are probably negligible. 

3.4. vg = 2 

3.4.1. Rotational data 

The rotational spectrum of methyl cyanide in its vg - 2 
excited vibrational state was also studied extensively. As for 
U 8 = l, the lower frequency data up to 112 GHz 1361 14611471 
were in part superseded by new data taken in Cologne; see 
Fig. [8] for an overview of the 7 = 2-1 transition. The remain¬ 
ing data were used in the final fit with reported uncertainties. 
The J" = 19 data from Ref. (37| had reported uncertainties 
of 50-150 kHz, though not specified. They fit on average to 
within 100 kHz and were retained in the final fit with assigned 
uncertainties of 100 kHz. The laser-sideband data |[38l were not 
retained in the fit for the same reasons as in the case of vg = 1 . 

Interactions between ug — 1 and 2 were described in subsec¬ 
tion [TT] 

The energy difference between vg = 2° and tq = 1 is 
203.54 cm -1 , rather large for a resonant Fermi interaction to 
occur. A AK = +1, Al = +2 interaction between K = 6 of 



Frequency / MHz 


Figure 8 : The rotational spectrum of CH 3 CN in the region of vs = 2, J = 
2-1 with largely resolved 14 N hyperfine splitting. All but one feature (near 
37010.5 MHz) used in the fit are indicated by assignment labels centered on the 
respective line. The first number gives the lower state F value. Primed F values 
indicate the weaker A F = 0 transitions; unprimed have A F = 1. Numbers after 
the comma signal the K value and, as superscripts, the / value. An ampersand 
was used to indicate overlapping HFS components. 



Frequency / MHz 


Figure 9: The rotational spectrum of CH 3 CN in the region of J = 75-74 
displaying the perturbation of i >8 = 2~ 2 , K = 6 by 1/4 = 1, K = 5. The K values 
are given for vs = 2 -2 ; the perturbation shifts the K = 6 transition lower by 
350 MHz. The clipped lines marked with a black pound sign are (from left to 
right) K = 4, 3, 2 and 1 transitions of vs = 1 +I . 


8 






































vg = 2 2 and K - 5 of u 4 = 1 was considered in the analysis of 
the V 4 band 1451 . We determined a resonance to occur at J = 73. 
Our closest ug = 2 transition is 7 = 75 - 74 which is calculated 
to be shifted by -350 MHz, see Fig.[9] A AK = ±3 (K = 8/5), 
Al = 0 interaction, resonant at 7 = 57, was identified here for 
the first time. Our closest vg = 2 transitions are J - 56 - 55 
and / = 59 - 58 with calculated shifts of about 20 MHz and 
25 MHz, respectively. 

Even though the energy difference between vg = 2 2 and v 1 = 
1 is 302.71 cm -1 , much larger than the difference between ug = 
2° and 14 = 1, the Coriolis interaction within the l components 
of ug = 2 2 and vj = 1 lead to a AK — 0, A/ = 3 Fermi resonance 
at K = 13 of ug = 2~ 2 and V 7 = 1 +1 , which caused perturbations 
in the vj band 1451 . We were able to assign transitions up to 
K = 11 in ug = 2 2 . Effects of this Fermi resonance were 
clearly present, but they were not sufficient to separate them 
reliably from the stronger Fermi resonance of ug = 2 2 with 
vg = 3 +1 occuring at K = 12. 

A AK = ±2, A1 = + 1 resonant interaction between K - 9 of 
ug = 2 ~ 2 and K = 1 of vj — 1~* near / = 88 shows no effects in 
the spectrum up to J - 65. At higher J , the intensities of these 
high energy lines are just too small to be identified reliably in 
our present spectral recordings. An identical type of interaction 
may occur at K — 14 and 15 of vg = 2° and K = 12 and 13 of 
U 7 = 1 +1 , however, we have no assignments yet for such high 
energy transitions. 

Fermi interactions between ug — 1 and 2 may also occur be¬ 
tween vg - 2 and 3. There are Aus = ±1, AK = 0, A 1 - +3 
Fermi resonances between i >8 = 2 2 and vg = 3 +I at K = 13 and 
between vg = 2° and ug = 3 +3 at K — 14/15; a resonant inter¬ 
action may occur for the latter at K = 15, currently predicted at 
J = 71/72. In the case of vg = 2“ 2 , our assignments extend to 
K = 11 and to K = 13 for vg = 2°. These levels are perturbed 
sufficiently, up to more than 20 MHz and almost 30 MHz, re¬ 
spectively, to evaluate the corresponding Fermi parameters. 

Ayg = ±1, AK = +2, Al — +1 interactions are expected to 
occur between vg = 2 and 3 at K = 12, l = -2 and K = 10, 
l = -1; at K = 14,1-0 and K = 12, / = +1; and at K — 16, 
/ = +2 and K = 14, Z = +3. One tentative assignment exists 
for the latter interaction. In addition, we have 6 assignments 
extending to as high as K = 19 of vg = 2 +2 , but several have 
residuals slightly larger than the estimated uncertainties, which 
may be explained by unregnonized overlap with other transi¬ 
tions, by perturbations of one or both rotational levels, or by 
underestimation of the uncertainties. 

We should point out that there may be interactions even 
higher in AK or A/, but we do not have evidence or data for 
these at present. No attempts were made to locate interations of 
vg — 2 with vibrational states higher than ug = 3. Such inter¬ 
actions may be non-negligible, but most likely at K levels with 
few or no assignments in our present data set. 

3.4.2. The 2vg band 

Transition frequencies were determined entirely from our 
new IR spectrum since Ref. Ii42l did not report their experi¬ 
mental values. 


The strongly allowed transitions in the 2vg IR band are those 
with AK = Al = 0. The K quantum numbers reach 12, higher 
by 2 than in the previous study l42l . while the maximum value 
of J is essentially identical, 66 versus 65. As in the previous 
study, we were unable to make any assignments in the perpen¬ 
dicular Al = ±2, AK = +1 component, which is only weakly 
allowed, but we were able to make numerous assignments to 
levels in the / = ±2 component which gain intensity through 
AK = Al = 2 mixing. This mixing is particularly strong for 
K — 4 of vg =2° and K = 2 of ug = 2 ~ 2 with the former being 
lower in energy at low-/ until the energy ordering changes at 
J = 40/41. As a consequence, many "RgU) and "P 4 (J) IR 
transitions were asssigned covering most of the J from 8 to 
62. It is interesting to note that the mixing of these K levels 
changes only slightly between adjacent J levels such that ro¬ 
tational transitions between these mixed K levels are too weak 
to be observed with the exception of those at the change of the 
energy ordering. As also shown in Table [2] the IR assignments 
to us = 2 ~ 2 cover mostly K — 1 to 3 with several K = 4 as¬ 
signments and few higher ones, those to vg = 2 +2 cover mostly 
K = 0 to 2; a few sparse higher-k' assignments are viewed as 
tentative and were weighted out in the final fit. 

Isolated lines with sufficiently good signal-to-noise ratio 
were given uncertainties of 0.0002 cm -1 , as in the case of vg. 
Noisier lines or lines in more crowded regions, such as in the 
(7-branch, were given larger uncertainties of 0.0003 cm -1 to 
0.0010 cnT 1 . 

We were also able to make assignments in the 3vs - V 8 hot 
band and in the V 7 - V 8 difference band, as was done earlier 
j42l . The latter band gains intensity because of an anharmonic 
resonance between 117 = 1 and vg = 3 sa. We will publish our 
detailed account on these states separately, fudging from the 
observation of transitions in 2v° with K = 12, we may be even 
able to make assignments in 4vg - 2vg. 

3.4.3. The 2vg - vs hot band 

Because of the strength of the 2vg band, little can be gained 
from the analysis of the 2v° - v x hot band. On the other hand, 
IR transitions to 2vg are limited to those which gain intensity 
through <722 interaction, see previous subsubsection. These tran¬ 
sitions belong to the / = -2 component for the most part. Ac¬ 
cording to Koivusaari et al. Il42ll . transitions of 2v 2 - vg are fre¬ 
quently overlapping each other or are blended with cold band 
transitions. Nevertheless, they made 52 r R -branch assignments 
with 3 < K” < 7. 50 of these assignments were made in the 
l = +2 component of 2vg - vs of which 6 had zero weight. 
Since no IR transition frequencies were published in Ref. Il42l . 
we searched our vs spectra for suitable hot band transitions. Un¬ 
fortunately, the congestion of the spectral lines and the quality 
of the spectra were such that only few lines could be assigned 
with certainty. 

3.5. Remarks on higher vibrational states 

The vibrational state ug - 2 interacts with tq = 1, with irj = 
1, and with ug = 3. Therefore, we need to take into account 
previously determined parameters for these states m. which 
involve interactions among these three vibrational states. 
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Table 3: Spectroscopic parameters or differences A thereof a,b (cm [ , MHz) c of methyl cyanide within vibrational states vs < 2. 


Parameter X 

« = () 

v% = 1 

d 8 = 2 U 

vs = 2 2 

£vib c 

0.0 

365.024 365(9) 

716.75042(13) 

739.148 225 (56) 

(A)* (A - B) 

148 900.103 ( 66 ) 

-115.930(26) 

-187.404(18) 

-259.956(122) 

(A)* B 

9198.899167(11) 

27.530277 (49) 

54.057 316(111) 

54.502729(70) 

(A)* D k x 10 3 

2 830.6(18) 

-11.46 (48) 

-20.19(34) 

-7.45(164) 

(A)* D jk x 10 3 

177.407 87(25) 

0.987 48 (60) 

1.675 5(25) 

1.808 8(15) 

(A)* Djx 10 6 

3 807.576(8) 

95.599(17) 

216.319(37) 

189.162(31) 

(A)* H k x 10 6 

164.6 ( 66 ) 

14.9 (22) 



(A)* H kj x 10 6 

6.062 0(14) 

0.034 1 (22) 

0.150 3(239) 

0.0249(55) 

(A)* H jk x 10 9 

1025.69(15) 

2.59(6) 

17.71(34) 

-0.37(21) 

(A)* Hj x 10 12 

-237.4(21) 

315.3 (30) 

200.7 (63) 

627.8 (59) 

(A)* L kk] x 10 12 

-444.3 (25) 




(A)* L jk x 10 12 

-52.75 (51) 




(A)* L jjk x 10 12 

-7.901 (32) 




(A f Ljx 10 15 

-3.10(17) 

-2.64 (20)^ 

-5.28(40/ 

-5.28(40/ 

(A)* P JK x 10 18 

552. ( 68 ) 




(A)* Pjjk x 10 18 

55.3(22) 




(A)* eQq 

-4.223 08(107) 

-0.038 7(19/ 

-0.077 4(38/ 

-0.077 4(38/ 

eQqrj 


0.1519(113) 



(A)* C bh X 10 3 

1.845 (90) 




(A)* ( C aa - C bb ) x 10 3 

-1.15(31) 




H 


138 656.195 (73) 


138656.042(102) 

riK 


10.3329(72) 


10.405 1 (98) 

BJ 


0.390469(7) 


0.394512(6) 

IKK X 10 6 


-834. (41/ 


-834. (41/ 

qjK x 10 s 


-34.057 (64/ 


-34.057 (64/ 

riJJ x 10 6 


-2.359 5(24/ 


-2.359 5(24/ 

qjKK x 10 9 


2.59(17/ 


2.59(17/ 

qjjK X 10 9 


0.508 5 (63/ 


0.508 5 (63/ 

q 


17.798438 (23) 

e 

17.729 857 (138) c 

q K X 10 3 


-2.664 5(111) 

e 

-2.615 3(89)'' 

qj x 10 6 


-63.842(14) 

e 

- 68.668 (31) e 

qjK X 10 9 


93.19(53/ 

e 

93.19(53/'" 

qjj X 10 12 


311.5(15) 

e 

191.9 (26)" 


a Numbers in pai*entheses are one standard deviation in units of the least significant figures. Empty entries indicate parameter not applicable or not used in the fit. 
See section [4] for sign and value considerations. 

b Parameter X given for u = 0; AX = Xi - Xo, with i representing an excited vibrational state. 
c All parameters given in units of megahertz, except for E’vib, which is given in units of inverse centimeters. 
d Ratio kept fixed in the fit for respective X or AX. 

e The parameter q and its distortion corrections connect levels with A K = A/ = 2, see subsection 13.11 and section[4] 
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There is a A K — Al — ±1 Coriolis interaction between 14 = 1 
and U 7 = 1 which will be resonant at K = 25 or slightly higher 
forv 7 = l -1 , too high to be observed easily. Another interaction 
occurs between K = 6/7 of 14 = 1 and K = 4/5 of 117 = 1 +1 . 

An anharmonic resonance occurs at K = 7/8 of 1)4 = 1 and 
Vg = 3 +3 . Another, stronger anharmonic resonance occurs at 
K = 7/8 of u 7 = 1 +1 and tig = 3 +1 . Inspection of preliminary 
energy levels involving these states suggest a plethora of further 
interactions among these states. Such interactions are also pos¬ 
sible with states of the next poly ad consisting of 1)4 = vg = 1 , 
14 = 1, U 7 = vg = 1, V(, - 1, and vg = 4 l69l . 

4. Determination of spectroscopic parameters 

Pickett’s spcat and spfit programs El were used for pre¬ 
diction of the CH tCN spectra and for fitting of the measured 
data. The programs were intended to be rather general, thus 
being able to fit asymmetric top rotors with spin and vibration- 
rotation interaction. They have evolved considerably with time 
because many aspects were not available initially GO ED, 
in particular special considerations for symmetric or linear 
molecules or for higher symmetry cases. Spectra of symmet¬ 
ric top rotors with excitation of one ElEl or two quanta of 
a degenerate vibrational mode E were fit successfully with 
spfit. The /-doubling is implemented as / = 0, ± 1 mod 3, thus 
/ = ±2 are identified with / = +1. Hence, G„, the A K = 0 Cori¬ 
olis interaction parameter between the /-components of 14 = 1 1 
and vg = 2? equals -2 AC, and +4,4 if, respectively. In fact, the 
Coriolis term enters the Hamiltonian as -2 (A£) v kl, see, e.g., 
Eq. (1) in Refs. If33l and ll42ll . The J- and /^-dependent dis¬ 
tortion parameters to G a are usually defined without the factor 
of -2, hence G J a equals r\j and -2r\j for vg - 1 and ug = 2, 
respectively, and so on. Similarly, e.g., G a and F ah , are the co¬ 
efficients of iJ a and (J a Jb + JhJa)/2, respectively. When they 
act as Coriolis parameters between two interacting states, the 
values from spfit often equal twice the values from other fit¬ 
ting programs. This is likely a consequence of the definition 
of G, which is approximately equal to 2//,/',. F designates a 
Fermi (more generally anharmonic) interaction parameter. J 
or K as sub- or superscripts represent, as usual, corresponding 
distortion corrections. Fi, Gi,, and F%j are the coefficients of 
operators which connect levels with AX higher by 2 than F, Gi, 
and Fjj, respectively, and with A/ modified accordingly to meet 
the symmetry condition A K - Al = 0 mod 3. I.e. Xi is defined 
for a parameter X as coefficient of {X, J + 2 L + 2 + J 2 k 2 ), with 
{a, b) = (ab + ba)/2 being the anticommutator. Interaction pa¬ 
rameters between different vibrations are indicated by numbers 
in parentheses separated by a comma; the degree of excitation 
and the / quantum numbers are also given if necessary. The 
most important interaction parameter with AK > 0 within a de¬ 
generate vibrational state of a strongly prolate symmetric rotor 
is q, also known as qi_i because it connects levels differing in 
AK — Al = 2. It can have large effects in particular in low-lying 
bending modes, in which the magnitude of C, often is close to the 
limiting case of 1 . 0 , which will cause ample near-degeneracies. 
Scaling factors for q depending on the /-components and the 
vibrational excitation are not implemented in spcat and spfit. 


The values are q/2 and ql V2 for Vg = 1 and 2, respectively. 
We note that q and its distortion corrections are sometimes de¬ 
fined negatively. This applies to the studies on methyl cyanide 
l33l l42l l45l l69l and propyne l50l [571 [52 1. Therefore, we in¬ 
verted the signs and added footnotes to the appropriate tables. 
We also note an additional factor of 2 in the definitions of q K , 
q JK , etc. in the articles on propyne, causing the determined 
numerical values to be smaller than in spfit. However, none of 
these parameters were used in the IR studies of methyl cyanide 
relevant to our investigation l33l [42;-45, [691. 

We determined rotational, centrifugal, and hyperfine struc¬ 
ture (HFS) parameters of the ground state as common for all 
vibrational states. As most of the data, in particular the IR data, 
were not affected by HFS, all states were defined twice, with 
and without HFS. Vibrational changes AX = X, - Xq, where 
X represents a parameter and X; the parameter in an excited 
vibrational state. This is very similar to several previous stud¬ 
ies on ch 3 cn mmmm and rather convenient because 
vibrational corrections AX are usually small with respect to X 
especially for lower order parameters X. In addition, this offers 
the opportunity to constrain vibrational corrections to vg = 2 to 
twice those of tig = 1, thus reducing the amount of independent 
spectroscopic parameters further. 

The degenerate vibrational states required AC, and q to be 
used in the fit with centrifugal distortion corrections as far as 
needed. Several of the distortion corrections between 14 = 1 
and 2 were constrained to be identical. Fitting of our new low 
frequency data of ug = 1 required eQqq to be used in the fit, 
the quadrupole parameter equivalent to q. We assumed that the 
value determined by spfit is +eQqi 7/4 

Starting values for the parameters of v = 0, vg = 1, and vg = 
2, were taken from Refs. EZHSJBZl. Parameters describing 
the triad around 1000 cm -1 , 04 = 1 , 07 = 1, and vg = 3, were 
taken from Ref. B31 and were kept fixed with the exception of 
one rotational interaction parameter between vg — 2 and 1)4 = 1 . 

Parameters were introduced newly to the fit, or a constraint 
was released, if the rms error (as the measure of the quality of 
the fit) was reduced by an amount deemed significant and if the 
new parameter was useful in the Hamiltonian and if the result¬ 
ing value(s) appeared to be reasonable. The latter aspect was 
a difficult one because vibrational corrections to some higher 
order parameters were of similar magnitudes as the parameters 
themselves and because some vibrational changes to vg = 2° 
or vg = 2? were rather different from twice the corresponding 
vibrational change to vg — 1. In most instances, care was taken 
to search for the largest reduction of the rms error. In almost 
all instances, this led to parameters determined with high sig¬ 
nificance, i.e. their values were more than five times larger in 
magnitude than their uncertainties. In few cases, parameters not 
determined with high significance were kept in the fit because 
the uncertainties were of similar magnitudes as those for other 
vibrational states. 

Concerning the interactions between vg = 1 and 2, we de¬ 
termined J distortion corrections Fj and Fjj to the main Fermi 
term F, whereas F K was not determined with significance. We 
estimated its value by assuming Fk/Fj * A/B. Once both F 2 
parameters were well-determined, their ratio was essentially 2 
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Table 4 : Spectroscopic parameters or differences A thereof a - b (cm 1 . MHz) c of methyl cyanide within vibrational states 04 = 1 , 1)7 = 1, and vg = 3 taken from 
Ref. ED. 


Parameter X 

1/4 = 1 

VJ = 1 

us = 3‘ 

us = 3 3 

£vib c 

920.29003 

1041.85471 

1077.786 3 

1 122.15 

A(A - B) 

-166.205 

889.440 

-302.030 

-347.76 

A B 

-46.148 22 

-5.73413 

80.294 85 

81.289 

A D k X 10 3 

-31.63 

149.333 

-22.8 rf 

-22.8 rf 

A D jk X 10 3 

7.165 

0.973 1 

3.504 &' 

3.504 6 rf 

AD j X 10 6 

-5.702 05 

14.001 

462. 

462. rf 

A H k X 10 6 

A H kj X 10 6 

A H jk x 10 9 

H 

-22.24 

-13.895 

315.68 

-148.1 

66 663.668 

138 665.87 

138 527.8 

VK 


7.224 8 

11.013 rf 

11.013 rf 

VJ 


0.078 153 

0.401 54^ 

0.401 54 d 

nKK X 10 6 
qjK X 10 6 

hjj x 10 6 


64.8 

-50.34 

2.385 

7.285 rf 

7.285 rf 

<1 


4.763 4 

17.683^ 

17.683 rf,c 

qj X 10 s 


-10.85 

-75.79 rf,e 

-75.79 rf ’ c 


a Empty entries indicate parameter not applicable or not used in the fit. Signs of q and qj altered, see section]?] 
b Parameter difference A(X) = X, - Xo, given for rotational and distortion parameters (7 represents an excited vibrational state). 
c All parameters given in units of megahertz, except for /T vl h, which is given in units of inverse centimeters. 
d Ratio kept fixed in the fit for parameter X or AX. 

e The parameter q and its distortion corrections connect levels with A K = A/ = 2, see subsection|3.1|and section]?] 


within the uncertainties. Therefore, the ratio was constrained 
to reduce the number of independent parameters. No distortion 
corrections could be determined for the Fermi terms associated 
with the interactions between v% — 2 and 3, so the J and K 
distortion corrections were fixed to appropriate multiples of the 
values determined for the v% = 1/2 resonance. Trial fits with Fi 
yielded values or constraints which were in some cases much 
smaller than those in the ug = 1/2 resonances with essentially 
no effect on the quality on the fit. These parameters were omit¬ 
ted from the final fits. 

The parameter /44 (for v% — 2) had little effects on the rms er¬ 
ror in most trial fits and was not determined with significance. 
It seemed to be determined in some trial fits, but its magni¬ 
tude appeared to be too large with respect to that of i>io = 2 
of propyne USD. Moreover, changes in some of the remaining 
spectroscopic parameters were deemed too large in these cases. 
Therefore, / 44 was omitted from the final fit. 

5. Discussion 

Our present ground state spectroscopic parameters are very 
close to those of our previous study ED; differences outside 
three times the experimental uncertainties occur only in Dj,Hj, 
and Lj. The experimental data associated with the v = 0/ug = 1 
perturbation help to constrain the K level structure in the ground 
vibrational state, additional contributions may come from the 
data involving the i>g = l/vg = 2 perturbations. As a conse¬ 
quence, we obtained a precisely determined value for Hk as 
(164.6 ± 6.6) Hz. The agreement with (156 ± 72) Hz from 
an earlier study l32l is good, though the large uncertainty in 
the earlier study limits the value of the agreement somewhat. 


Table 5: Interaction parameters' 7 (MHz) between vibrational states of methyl 
cyanide. 


Parameter 

Value 

F 2 (0,8*)x 10 3 

-70.897 (27) 

F(8 ±1 ,8 2>:p2 ) 

53 157.7(33) 

F k ( 8 ± 1 , 8 2 ’* 2 ) 

- 6 . b ’ c 

Fj( 8 ±1 ,8 2 - :f2 )x 10 3 

-369.89 (44)* 

Fyy(8 ±1 ,8 2:F2 )x 10 6 

1.681 (87)* 

F 2 (B ±1 ,8 2 '°)x 10 3 

-65.491(24/ 

F 2 (B ±1 ,8 2 - ±2 )x 10 3 

-130.982(48/ 

FaA 8 2 - ±2 ,4‘) [2VV488] 

8.7362(21) 

F 2ac (S 2 -°,4‘)xl0 6 

-7.98(16) 

F( 8 ± 2 , 7 ^) [W788] 

45 170.8^ 

F(8 2,±2 , 8 3 ’ t1 ) 

77 208.(93) 

F(8 2 ’ 0 ,8 3 - 3 ) 

91509.(131) 

G*(4.7) [ 2 W 47 ] 

909. e 

FbA 4,7) [2w 47 ] 

-1.84 f 

F(4,8 3 ' ±3 )) [W4888] 

11 430. e 

F’(7 ±1 ,8 3 ±1 ) [3^7888] 

50 129.2 C 

GairW-^VW^g} 

-2239.T 


a Alternative designations from Ref. ED given in brackets. Numbers in 
parentheses after the interaction parameter designate the vibrational states 
separated by a comma, see also section [4] Numbers in parentheses after the 
values are one standard deviation in units of the least significant figures. 
b J and K distortion corrections to / r (8 2,±2 , 8 3,:f1 ) and F(8 2,0 ,8 3,3 ) kept fixed 
to V2 and V3, respectively, times the corresponding F(8 ±1 ,8 2 value; see, 
e.g., Ref. EQ. 

c Estimated assuming Fr/Fj « A/B, see section[4] 
d Ratio constrained, see section [ 4 ] 
e Kept fixed to values from Ref. ED 


12 










Table 6: Spectroscopic parameters or differences A thereof a,b (cm 1 , MHz) c of methyl cyanide involving v& = 1 and 2 in comparison to previous values. 


Parameter X 

y 8 

= 1 

y 8 


vs 

= 2 1 

Present 

Kef. (33] 

Present 

Ref. 1421 

Present 

RefTl42l 

£vib c 

365.024365 (9) 

365.015965 (12) 

716.75042(13) 

716.70441 (1) 

739.14823 (6) 

739.12907 (3) 

AA 

-88.400 (26) 

-90.079 (36) 

-133.347 (18) 

-140.069 (39) 

-205.453 (122) 

-210.826 (96) 

A B 

27.53028 (5) 

27.53255 (33) 

54.05732(11) 

54.07405 (54) 

54.50273 (7) 

54.50740 (90) 

A D k X 10 3 

-11.46 (48 ) 

-18.11 (81) 

-20.2 (3) 

-12.3 (10) 

-7.5 (16) 

20.1 (42) 

A D jk X 10 3 

0.9875 (6) 

0.8739 (129) 

1.6755 (25) 

1.3316(69) 

1.809(2) 

1.790(39) 

A Dj X 10 6 

95.599(17) 

93.436 (69) 

216.32(4) 

218.41 (18) 

189.16(3) 

182.93 (24) 

A H k X 10 6 

14.9 (22) 

-215.0 (51) 


-384.3 (72) 



A H kj x 10 6 

0.034 (2) 

0.908(114) 

0.150 (24) 


0.025 (6) 


A H jk x 10 9 

2.59 (6) 

35.74(171) 

17.71 (34) 


-0.37 (21) 


A Hj X 10 12 

315.3 (30) 


200.7 (63) 


627.8 (59) 


A Lj X 10 15 

-2.64 (20/ 


-5.28 (40/ 


-5.28 (40/ 


A(x 1CT 3 

138.65620 (7) 

138.64715 (7) 



138.65604(10) 

138.63733 (16) 

VK 

10.333 (7) 

10.448 (93) 



10.405 (10) 

10.448 (9) 

VJ 

0.390469 (7) 

0.389679 (66) 



0.39451 (1) 

0.39420(18) 

VKK X 10 6 

-834. (41/ 

-453. (174) 



-834. (41/ 


t]j K x 10 6 

-34.06 (6/ 

-40.77 (90) 



-34.06 (6/ 

-46.8 (30) 

qjj x 10 6 

-2.3595 (24/ 

-2.6750(111) 



-2.3595 (24/ 

-2.845 (45) 

>1KKK X 10 6 


15.23 (90) 





r/JKK X 10 9 

2.59 (17/ 




2.59(17/ 


' lJJK X 10 9 

0.509 (6/ 




0.509 (6/ 


1 

17.79844 (2) 

17.79776(11) 

e 

e 

17.72986 (M)' 

17.70068 (51/ 

q K X 10 3 

-2.6645 (111) 


e 

e 

-2.6153 (89) c 


qj X 10 6 

-63.84(1) 

-63.80(16) 

e 

e 

-68.67 (3) e 

-66.84(17/ 

qjK X 10 9 

93.19 (53/ 


e 

e 

93.19 (53/- e 


qjj X 10 12 

312.(2) 

324. (57) 

e 

e 

192. (3 y 



a Numbers in parentheses are one standard deviation in units of the least significant figures. Empty entries indicate parameter not applicable or not used in the fit. 
See section [4] for sign and value considerations. 

b Parameter difference A(X) = X; - Xo given for rotational and distortion parameters (i represents an excited vibrational state). Signs were adjusted to present 
definition for data from Refs. >3311421 . In all instances, previous values refer to the combined fits of IR and rotational data. 
c All parameters given in units of megahertz, except for E v fo, which is given in units of inverse centimeters. 
d Ratio kept fixed in the fit for respective X or AX. 

e The parameter q and its distortion corrections connect levels with A K = A/ = 2, see subsection 13.1 | and section[4] 
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Our previous estimate of 51 Hz, derived by scaling D K with the 
Dk/A ratio as an order-of-magnitude estimate, was apparently 
too small. 

As can be seen in Table© the extensive new rotational data 
pertaining to v& - 1 and 2 allowed us to determine spectro¬ 
scopic parameters or changes thereof up to higher order than 
previously MM- The uncertainties of many lower order pa¬ 
rameters were improved in addition. 

The previous and present spectroscopic parameters agree rea¬ 
sonably at lower order and particularly for those with depen¬ 
dence on J. The agreement deteriorates for higher order pa¬ 
rameters, especially those dependent on K. Consideration of 
the resonances between us = 1 and 2 as well as between v% = 2 
and 3 affected the purely K -dependent parameters significantly, 
very similar to the case of propyne, where the Fermi interac¬ 
tions involving the equivalent vio band occur at slightly lower 
K J50)(5TJ. As in the case of propyne, we did not require tikkk 
to be included in the fit. Comparison of the various q param¬ 
eters from our study as well as those from Ref. 1331 suggests 
that the magnitude of ijkkk reported in that study is at least a 
factor of 100 too large. Similarly, the previous A Hk value of 
i;8 = 1 in CH 3 CN was -215 ±5 Hz |33) . whereas our value is 
14.9 + 2.2 Hz; in addition, a value of -384 + 7 Hz was reported 
for us = 2°, in Ref. 1331 . whereas we obtained a value similar to 
that in vg = 1 in trial fits. However, the uncertainties were only 
slightly smaller than the values, so AHk(v& = 2 °) was omitted 
in the final fit. 

The vibrational changes in distortion parameters from v = 0 
to us = 1 in Tableware often small, though not in all instances. 
The rather large changes in Hj and Lj are common in these 
types of molecules. In fact, the changes of Hj upon excitation 
of each of the two low-lying bending modes are even larger 
in magnitude in propyne EE). Excitation of the CCN bending 
mode leads also to large changes in H of HC3N 11771 . DC3N 
1781 . or other, similar molecules. 

The changes of rotational and centrifugal distortion parame¬ 
ters from ground state values in v% = 2 ° and us = 2 1 are roughly 
twice those in u$ = 1 , though pronounced deviations occur al¬ 
ready in AA (or A(A - B)). The similarity to related values 
in propyne fl50l I5T1 and methyl isocyanide 17511761 may indi¬ 
cate that the low energies of the vibrations lead to substantial 
changes in A A with different / quantum number. Changes upon 
excitation in i >8 are closer to the ratio of 2 in A B, but devia¬ 
tions increase, in particular for the changes in the AH values. 
These deviations may be caused by the incomplete analysis of 
the us = 2/3 or other resonances. 

The A A values are essentially identical between V& = 1 and 
i )8 = 2?. The changes in i]k and rjj are rather small, and higher 
order parameters were constrained to the same values. The q 
values are rather similar also, whereas the qjj values differ con¬ 
siderably, possibly again a consequence of a not or not fully 
accounted interaction. 

With respect to the perturbation of the ground vibrational 
state by t )8 = 1, Simeckova et al. ED calculated perturba¬ 
tions of about 45 MHz and 36 MHz for the two ground state 
transitions involving / = 43 and K — 14, this is about 20 MHz 
larger than actually observed. The difference can be traced al¬ 


most entirely to the energy differences between the two perturb¬ 
ing levels at / = 43, which was estimated to be 0.013 cm -1 
ED whereas our calculations yield 0.025 cm 1 . Our value 
of F flr (8 2,±2 ,4 l )/2, 4.368 ± 0.001 MHz agrees quite well with 
|w 488 l = 4.101 ± 0.042 MHz from the combined analysis of IR 
and rotational data or with 4.170 + 0.042 MHz from an analysis 
of IR data only ED- We should point out that the sign cannot 
be determined from the fit of the line frequencies, but it may 
be determinable from relative intensities in the rotational or IR 
spectra. The value given in Ref |45| was actually negative; we 
will try to determine the sign in our subsequent study which 
will include 1)4 = 1 data. Other interaction parameters between 
different vibrational states were either determined for the first 
time here or were not released in our fits. 

Selected low order spectroscopic parameters of CH 3 CN 
have been compared in Table |7] with corresponding values of 
CH 3 CCH and CH 3 NC. The data of CH 3 CN and of CH 3 CCH 
are comparatively similar, while those of CH 3 NC with its much 
lower bending mode are substantially different. The Fermi pa¬ 
rameters between the first and the second bending state reflect 
these changes, and the values appear to scale approximately 
with q\,. 

F( 8 " ±2 , 8 3,:f1 ) and F( 8 2 , 0 , 8 3 - 3 ) should ideally differ from 
F( 8 ±l , 8 2,:f2 ) by factors of V 2 and V3, respectively, see, e.g.. 
Ref. J5D- The experimental values of 77208 ± 93 MHz and 
91509 + 131 MHz compare quite well with the expected values 
of 75192 MHz and 92091 MHz, respectively, in particular con¬ 
sidering that there are no i >8 = 3 transition frequencies in the 
fit, and that we used spectroscopic parameters from a previous 
study for this state fl45l . 

The lack of assignments in 2v8 - vs may limit the accuracies 
of the predicted line positions in the 2 vg 2 - Vg 1 sub-band at 
higher K. However, our data on the perturbations between ug = 
1 and 2 may in part compensate these limitations. 

6. Expanded database for CH3CN line parameters 

Compilations of CH 3 CN molecular line parameters are 
greatly expanded by including new or updated predictions of 
spectra within or between the three states of the present in¬ 
vestigation. Predictions of rotational and infrared spectra for 
astrophysical applications will be available in the catalog sec- 
tiorjjof the CDMS Q)E) and that of jpifltza. The parameter, 
line, fit and other auxiliary files are available online as supple¬ 
ments with the publication and in the Cologne Spectroscopy 
Data sectiorj^Jof the CDMS. For remote sensing of, e.g.. Titan, 
the new compilation will also be provided to HITRAN lt64l and 
GEISA lf80l with ^-pressure and self-broadening information 
included. For this, the extensive pressure broadening CH 3 CN 
measurements of P and R branch transitions of V 4 l 68 l were 
modeled by extending empirical expressions applied to ammo¬ 
nia 18111821 . The resulting polynomial has no physical meaning, 

4 http://www.astro.uni-koeln.de/cdms/catalog 

5 http://spec.jpl.nasa.gov/ftp/pub/catalog/catdir.html 

6 http://www.astro.uni-koeln.de/site/vorhersagen/daten/CH3CN/CH3CN/ 


14 



Table 7: Comparison of selected low order spectroscopic parameters or dif¬ 
ferences A thereof 1 (cm -1 , MHz ) 6 of methyl cyanide, methyl acetylene, and 
methyl isocyanide in their ground and lowest bending states b. 


Parameter X 

CH 3 CN c 

CHjCCH" 

CH 3 NC c 

Ao 

158 099. 

159139. 

157 308. 

B 0 

9198.9 

8 545.9 

10052.8 

Dk.o x 10 3 

2831. 

2939. 

2 567. 

Djk.o x 10 3 

177.4 

163.4 

227.5 

Dj ,0 x 10 6 

3 808. 

2908. 

4692. 

E b b 

365.0 

330.9 

267.3 

A 11 A b 

-88.3 

-61.9 

0.6 

A“S b 

27.53 

23.94 

39.07 

Mb 

138 656. 

141919. 

146 007. 

VK.b 

10.34 

10.94 

10.13 

VJ .b x 10 3 

390.5 

334.8 

587.6 

lb 

17.80 

16.79 

6.94 

Ffb* 1 , b 2 ’ T2 ) 

53 169. 

51745. 

23 808. 


0 Parameter X given for v = 0; A(X) = X, - Xo, with i representing an excited 
vibrational state. 

b All parameters given in units of megahertz, except for E b , which is given in 
units of inverse centimeters. 
c This work. 

cl Ref. 151 1 ; sign of q adjusted, see sectionjd] 
e Refs. 153117511761 : sign of q adjusted, see section[4] 


but it does provide reasonable estimates of Lorentz broadening 
coefficients for CH 3 CN transitions. 

6.1. Pressure-broadened half widths 

The measured N 2 - and self-broadened halfwidths (HWHM) 
for CH 3 CN in the V 4 band at 920 cm -1 vary smoothly with re¬ 
spect to their quantum numbers J and K ll 68 l . Patterns become 
more obvious when the widths are examined as a function of 
values related to the transition quantum numbers; these are de¬ 
fined as J m = J" + 1, K m = K" + 1 for R branch transitions, and 
J m = J", K m = K" for P and Q branch transitions. Examples 
are shown in Table |8) N 2 - and self-broadened half widths were 
derived for our CH 3 CN compilation by modeling the reported 
widths for the V 4 band |} 68 l according to the expression in Eq.[T| 
the polynomial coefficients are listed in Table [9] 

y(J m , K m ) = £ Cl U J in K L (1) 

hi 

The measured widths were weighted in the fit by the 
inverse of the square of the reported uncertainties 8681 . 
The impact of a few outliers (at very low and higher J m 
and K m ) were suppressed as an important side-effect. As 
shown in Fig. [TO] the calculated values using Eq. |T| (solid 
curves) reproduce the pattern of the measured data points 
within 5 % for most of the transitions. The full lists are 
given in Supplemental files 5 Width(N2) CH3CN v4.out and 
6 Widthf Self) CH3CN v4.out. Finally, we adopted this expres¬ 
sion to estimate N 2 - and self-broadened half widths for the three 
IR bands, vg, 2v s - vg, 2vg of the present study. Because the 
experimental V 4 lines, for which line widths were determined, 
have J m < 50, we limited the line-widths at 296 K to: 


Table 8 : Examples of measured 0 (obs.; with uncertainties unc. in percent in 
parentheses) and modeled (calc.; with percent difference in parentheses) V 4 self- 
broadened line widths (cm _l /atm) of methyl cyanide at 296 K with frequencies 
(cm -1 ) and quantum numbers. 


J'.K' 

J", K" 

Frequency 

Self-broadened width 
obs. (%unc.) calc. (%diff.) 

Jm 5 X m 

2, 1 

3, 1 

918.43291 

1.358(0.7) 

1.2588 (7.88) 

3,1 

3, 1 

2, 1 

922.10564 

1.307(1.0) 

1.2588 (3.83) 

3, 1 

7,6 

8,6 

915.09499 

1.071 (0.7) 

1.0798 (-0.82) 

8,6 

8,6 

7,6 

924.88187 

1.039 (0.9) 

1.0798 (-3.78) 

8,6 

18,9 

19,9 

907.67687 

1.722 (0.5) 

1.6957 (1.55) 

19, 9 

19,9 

18,9 

930.89430 

1.711 (0.7) 

1.6957 (0.90) 

19, 9 

22, 7 

23, 7 

905.12404 

1.773 (0.4) 

1.7497 (1.33) 

23, 7 

23,7 

22. 7 

933.24349 

1.821 (0.6) 

1.7497 (4.07) 

23,7 

28,0 

29,0 

901.25148 

1.583 (0.5) 

1.5807 (0.15) 

29, 0 

29,0 

28,0 

936.73087 

1.528 (0.5) 

1.5807 (-3.33) 

29,0 

40,4 

41,4 

892.54529 

0.896 (0.7) 

0.8737 (2.55) 

41,4 

41,4 

40,4 

942.65553 

0.823 (0.8) 

0.8737 (-5.81) 

41,4 


a The measured values are taken from Ref. (Si. 


Table 9: Coefficients to model the N 2 - and self-broadened widths of the V 4 band 
of CH 3 CN using Eq.|T] 


Coefficients 

n 2 

Self 

ooo 

0.171 

1.621 

a\Q x 10 3 

4.731 

-172.4 

tfoi x 10 3 

-9.962 

-82.99 

a 20 X 10 3 

-0.7152 

26.32 

a\\ x 10 3 

1.704 

18.89 

<202 x 10 3 

-1.233 

-14.42 

^30 x 10 3 

0.03155 

-1.28 

«21 X 10 3 

-0.1068 

-1.296 

a \2 x 10 3 

0.1361 

0.7237 

<203 x 10 3 

0.002544 

0.1243 

040 x 10 6 

-0.6094 

24.83 

<331 x 10 6 

2.643 

-34.66 

<322 X 10 6 

-3.849 

-10.5 

a\ 2 , x 10 6 

-1.819 

9.179 

<350 x 10 9 

4.388 

-171.6 

<341 x 10 9 

-22.42 

-319.8 

032 X 10 9 

33.41 

66.99 

023 X 10 9 

37.85 

-430.9 
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0.12 < y N , < 0.2 cm '/atm for N 2 -broadened half-widths, 
and 

0.2 < y se ]f < 2.0 cm - 1 /atm for self-broadened half-widths. 

We adopted one single value, 0.72, from the HITRAN2012 
database iRyfl for the temperature dependence exponent of the 
No-broadened Lorentz widths. This value is based on the mea¬ 
surements of the J - 12-11, K = 3 transition in the ground 
vibrational state lf83l . 


6.2. N 2 -pressure-induced frequency shifts 

Rinsland et al. (681 also reported N 2 -pressure-induced fre¬ 
quency shifts for CH 3 CN in its V 4 band, but no distinctive pat¬ 
tern is seen for pressure shifts with respect to J,„ and K m . How¬ 
ever, there was a significant correlation between measured line 
widths and frequency shifts, as presented in their Fig. 11 of 
Ref (68l . We use this observation in our present work to obtain 
frequency shift estimates for the transition of vg, 2 v s - vg, and 
2vg. First, we apply a scaled set of the measured N 2 -induced 
frequency shifts for V 4 transitions for all other CH 3 CN tran¬ 
sitions having the same set of quantum numbers, J", K ", J', 
and K'. However, not all transitions predicted in the present 
work have counter parts among the measured V 4 transitions; we 
assume for these cases that their N 2 -induced frequency shifts 
64 (J m , K m ) can be estimated from their pressure-broadened line 
widths using the width-shift relationship proposed by Rinsland 
et al. [68 j. Based on Fig. 11 of Ref. l l68ll . a simple expression 
was derived as 


<y 4 (/», K,n) = [yiJm,Km) n 2 - 0.15J/15 (2) 

One commonly assumes that the overall magnitudes of the 
pressure shifts increase with increasing vibrational band centers 
(84}. We estimated the pressure shifts 6 fJ m ,K m ) of v, from 
those of the V 4 band according to 


6.3.1. v 8 

Our simulations on a line-by-line basis indicate that the vg 
spectrum is represented reasonably well by a transition dipole 
moment of p( 8 g) = 0.043 D. 

We calculated partition function values by summing over all 
states up to / = 99, K — 25, and vg — 3. While the rotational 
contribution to the partition function at 296 K is well converged, 
the vibrational contribution is slightly truncated. We have esti¬ 
mated the truncation error classically, assuming a rigid rotor as 
well as a harmonic oscillator. This approximation in itself in¬ 
troduces some error, but we estimate it to be small with respect 
to the truncation error. The truncation leads to a corrections 
of 1.408%. The resulting intensity, integrated over all lines, 
is given in Table 10 The integrated intensity in the entire vg 
region, i.e. with all hot bands, was estimated by multiplying 
the calculated vg intensity with the vibrational partition factor 
of 1.5006. An additional increase of the integrated intensity by 
2.677 % was required because of the contributions of other iso¬ 
topic species. There is a convention in HITRAN that intensities 
refer to a sample in natural isotopic composition. The resulting 
intensity is also given in Table 


between two lower values (67l 


Our value is about half way 
and one higher one | 


6.3.2. 2vg - vg 

a/ 

We assumed transition dipole moments of /A 8 y) = 0.043 D 
for both / = 0 and / = +2 components as for the cold band. This 
approximation appears to be valid well within the estimated in¬ 


shows the calculated contributions of the vg and 2 vg - vg sub¬ 
branches. 

As can be seen in Table [TO] slightly more than 10% of the 
total intensity in the vg region of methyl cyanide are caused by 
IR transitions originating in states higher than t>g = 1. Finding 
unblended lines among these will be difficult because this was 
already not easy for 2 v\ - vg. 


tensity uncertainties of 4 %, see subsubsection [TO] Fig. m 


6j(J m , K m ) = 64 (J„„ K m ) x (E vib (0/£vib(4)) (3) 

with FVibO) being the vibrational energy of vg, 2 vg - vg, or 
2vg and £’ v ib(4) that of V 4 . 

The resulting ^-pressure shifts for the vg, 2vg - vg, and 2vg 
bands are presented in Fig. |T7] i. The width-shift relationships 
resemble those observed in the V 4 band ( 68 ]. The distribution 
with respect to K m , shown in Fig. [IT)), does not show a distinc¬ 
tive pattern, as was the case for Vi | 68 j. 


6.3. Intensities of v g, 2vg - vg, and 2vg 

We evaluated transition dipole moments for each of the three 
IR bands of the present study. There are some indications that 
additional correction terms may lead to a better intensity mod¬ 
elling, especially in the cases of vg and 2vg - vg. More so¬ 
phisticated intensity modelling is beyond the aim of the present 
study, in particular because the quality of the IR spectrum in the 
vg region is not sufficient for such modelling. 


6.3.3. 2vg 

The relative precisions of measured intensities were checked 
using spectral comparisons. For example, the 731-732 cm -1 
interval in Fig.@ has 181 retrieved features with intensities 
ranging from 1.3 x 10 6 to 7.3 x 10 3 cm 2 /atm at 293 K. Since 
the sample was a mixture, line intensities of the two calibration 
species were measured to obtain the partial pressures of OCS 
(3 %) and CO 2 (0.4 %). In addition, selected strong CH 3 CN 
lines of the V 4 band were retrieved and compared to values re¬ 
ported in Table 2 of Ref. ( 68 l ; this indicated that the retrieved 
intensities in the 2vg region should be increased by 4%; with 
these calibrations, the overall (relative) precision is thought to 
be 0.0002 cm -1 for line centers and 5 % for intensities. How¬ 
ever, worse precisions are expected for blended features sep¬ 
arated by less than 0.0025 cm -1 and for ~5100 absorptions 
weaker than 9 x 10 5 cm 2 /atm at 293 K. 

Our simulations indicate that the spectrum is represented 
well by a transition dipole moment of p( 8 g°) — 0.030 D. See 
subsection 2.3 for considerations concerning the experimental 
intensities and subsubsection 6.3.1 for a remark on the partition 
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Figure 10: Model fit of measured Ns- (left) and self-broadening (right) Lorentz width coefficients in cm '/atm at 296 K for the V 4 band ED- Many of the widths 
are well reproduced by Eq.|T] The solid lines in green represent calculated values while observed widths of P and R branch transitions are indicated by circles in 
red and hexagons in blue, respectively. In the lower parts, the J m and K m values are combined for the horizontal axes to reveal the variation patterns. 


(a) 0.015 


_ 0.01 

E 



-O.OIS foj^ig- f * friiMSr 

0.1 0.12 0.14 0.16 0.18 0.2 
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Figure 11: (a) ^-pressure induced frequency shifts derived for the vs region. The scattered values are adapted from the V 4 measurements (Ml while the well-aligned 
ones are based on Eqs. [2]and[3] (b) The derived N 2 -shifts of IR bands in the present study are presented for selected K m (see text for details). 


Table 10: Derived transition dipole moments /i (D) and integrated intensities (10 -19 cm/molecule) of IR bands of 12 CHpC 14 N calculated from the present Hamil¬ 
tonian model with respect to the transition dipole moments (Calc. (jj .)) and apparent intensities in the / components (Calc. (/)) together with the extrapolated and 
experimentally measured total integrated intensities from previous studies. 


IR band 

n 


Integrated intensity 



Calc. (//) 

Calc. (/) 

Ref. ( 66 J 

Ref. (67J 

Ref. 1851 

Ref. ( 86 ) 

I>8 = 1 — 0 

0.43 

1.483 

1.483 





V S = 2° - 1 

0.43 

0.251 

0.244 





I>8 = 2 2 - 1 

0.43 

0.257 

0.264 





vs region 


2.285° 



1.77(4) 

1.81 

2.79 (28) 

I>8 = 2° - 0 

0.30 

1.674 

1.585 





«8 = 2 2 - 0 

0.00 

0.000 

0.090 





2 v 8 region 


2.580° 


2.63 (18) 

2.50 (4) 

3.43 



a Estimated from the cold bands, ug = 1-0 and 2 - 0, by multiplication with the vibrational partition factor 1.5006 at 296 K. An increase of 2.677 % from other 
isotopic species was considered. 
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function and on the sample composition. The integrated inten¬ 
sity from our Hamiltonian model is 1.75 x 10 19 cm/molecule 
for the 2 v 8 band alone. Only a small portion of this intensity, 
~ 1.0 x ICE 21 and ~ 8.1 x ICE 21 cm/molecule is transfered 
through the <722 interaction to the vg = 2 +2 and vg = 2 2 compo¬ 
nents, respectively. 

About 23 % of the 2vg intensity at 296 K should arise from 
the 3 v 8 - vg and V 7 - vs bands. A preliminary analysis suggests 
that a fair fraction of these lines will not be blended by other 
lines of similar or higher intensity. Even some of the lines orig¬ 
inating in states higher than vg = 1 , which account for slightly 
more than 10 % of the band intensity, may be usable for spectro¬ 
scopic analyses. More work is needed to clarify these aspects. 

7. Conclusions 

Our extensive and accurate rotational data in vibrational 
states vg <2 resulted in a parameter set extended to higher order 
and to mostly improved uncertainties at lower order. More im¬ 
portantly, these data enabled the first treatments of resonances 
between ug = 1 and 2, as well as between vg — 2 and 3. These 
resonances cause considerable perturbations in the rotational 
spectrum, and their successful treatments affect the spectro¬ 
scopic parameters, in particular the /(-dependent ones and those 
at higher order. We have also performed the first in-depth anal¬ 
ysis of the resonant interaction between the ground vibrational 
state and vg — 1 . 

The analysis of a new 2vg spectrum was an additional im¬ 
portant contribution to the determination of spectroscopic pa¬ 
rameters. Only about half of the roughly 12000 lines in the 2vg 
spectrum are caused by the cold band. A large fraction of the 
unaccounted lines will be caused by the 3vs - vs hot band and 
the v 7 - vg difference band. Accounting for these will require 
in depth analysis of the v' 7 /3vg band system aided by rotational 
transitions. A still considerable fraction of lines are due at least 
in part to bands originating in states higher than vg — 1 . 

We have derived line broadening and shift parameters for the 
vs, 2 v 8 -V 8 , and 2 vs bands from previously determined V 4 values 
and have used this information to determine transition moments 
and intensities for the three IR bands. Individually measured 
transition intensities may lead to an improved intensity model. 

We have generated an expanded database for CH 3 CN line pa¬ 
rameters for transitions within each of the states vg < 2 , mainly 
for radio astronomical observations, in particular with ALMA 
and other interferometers. In addition, information is provided 
for transitions between these three states for remote sensing of, 
e.g., the atmosphere of Titan. 

The inclusion of data pertaining to 14 = 1 is almost complete. 
We intend to improve the data involving its interactions with 
other vibrational states. We want to study also the proposed 
strong, although distant, thus non-resonant Fermi interaction of 
14 = 1 with vg = 2 ||43| . In addition, we have already made 
assignments in the vq + vg - V 8 IR band as well as for rotational 
transitions in 04 — vg = 1, though these are often tentative. The 
spectra are perturbed, but, as it appears, only locally, such that a 
reasonable deperturbation may be possible by considering only 
spectroscopic parameters of the perturbing states. 



Figure 12: Contributions of each sub-band in the vs region of CH 3 CN showing 
their relative intensities for the cold band in the upper part and for the hot band 
in the lower part. Note: hot bands originating in states higher than v& = 1 
contribute to ~ 10 % to the total intensity. 
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Appendix A. Supplementary Material 

Supplementary data for this article are available on Sci- 
enceDirect (www.sciencedirect.com) and as part of the 
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data associated with this article can be found, in the online ver¬ 
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readme file are provided as well as the full tables of the model¬ 
ing ofN 2 - and self-broadening line widths of the V 4 band along 
with an explanation file. 
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